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Abstract 



In the context of one-dimensional diffusions, we present basic esti- 
mates (having the same lower and upper bounds with a factor of 4 only) 
»^ ' for four Poincare-type (or Hardy-type) inequalities. The derivation of 

, two estimates have been open problems for quite some time. The bounds 

provide exponentially ergodic or decay rates. We refine the bounds and 
illustrate them with typical examples. 

> 

. 1 Introduction 

An earher topic on which Louis Chen has studied is about the Poincare-type 
inequalities (see [B [2] , for instance) . We now use this good chance to introduc- 
tion in Section [2] some recent progress on the topic, especially for one-dimen- 
sional diffusions (elliptic operators). The basic estimates of exponentially er- 
godic (or decay) rate and the principal eigenvalue in different cases are pre- 
sented. Here the term "basic" means that upper and lower bounds are given 



^ • by an isoperimetric constant up to a factor four. As a consequence, the criteria 



for the positivity of the rate and the eigenvalue are obtained. The proof of the 
main result is sketched in Section [3l The materials given in Sections HI [Sj and 
Appendix are new. In particular, the basic estimates are refined in Section U] 
and the results are illustrated through examples in Section [5j The coincidence 
of the exponentially decay rate and the corresponding principal eigenvalue is 
proven in Appendix for a large class of symmetric Markov processes. 
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2 The main result and motivation 

2.1 Two types of exponential convergence 

Let us recall two types of exponential convergence often studied for Markov 
processes. Let Pt{x,-) be a transition probability on a measurable state space 
{E, S) with stationary distribution vr. Then the process is called exponentially 
ergodic if there exists a constant e > and a function c(x) such that 

\\Pt{x, •) - ^llvar ^ c{x)e-'\ t^O, xeE. (1) 

Denote by Smax be the maximal rate e. For convenience, in what follows, we 

allow £max = 

0. Next, let L'^{7t) be the real L2(7r)-space with inner product 
(•, •) and norm || • || respectively, and denote by {Pt}t^o the semigroup of the 
process. Then the process is called to have L'^ -exponential convergence if there 
exists some r] 0) such that 

\\Ptf - vr(/)|| < 11/ - 7r(/)||e-''*, t ^ 0, f G L^tt), (2) 

where vr(/) = /dvr. It is known that r/max is described by Ai: 

Ai = inf{(/, -Lf) : / e ^(L), n{f) = 0, ||/|| = 1}, (3) 

where L is the generator with domain S'{L) of the semigroup in L^{tt). Even 
though the topologies for these two types of exponential convergence are rather 
different, but we do have the following result. 

Theorem 2.1 ([3], |6|) For a reversible Markov process with symmetric measure 
TT, if with respect to vr, the transition probability has a density pt{x,y) having 

1/2 

the property that the diagonal elements ps{-, •) G Lj^^ (tt) for some s > 0, and a 
set of bounded functions with compact support is dense in L^(vr), then we have 

'^max — Ai. 

As an immediate consequence of the theorem, we obtain some criterion for 
Ai > in terms of the known criterion for emax > 0. In our recent study, we 
go to the opposite direction: estimating ffmax in terms of the spectral theory. 

We are also going to handle with the non-ergodic case in which ([2]) becomes 

KiPtff) < t ^ 0, / G L\f,), (4) 

where is the invariant measure of the process. Then r/max becomes 

Ao = inf { - /i(/L/) : / e f,{f) = l}, (5) 

where is a suitable core of the generator, the smooth functions with compact 
support for instance in the context of diffusions. However, the totally vari- 
ational norm in ([T]) may be meaningless unless the process being explosive. 
Instead of ([1]), we consider the following exponential convergence: 

Pt{x, K) ^ c{x, K) e"^*, t ^ 0, x£ E, K: compact, (6) 
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where for each compact c{-,K) is locahy ^-integrable. Under some mild 
condition, we still have ffmax = Aq- See Appendix for more details. 

2.2 Statement of the result 

We now turn to our main object: one-dimensional diffusions. The state space 
is E := {—M,N) {AI,N ^ oo). Consider an elliptic operator 

where a > on Then define a function C{x) as follows: 

r b 

C{x) = / -, xe E, 
Je a 

where 6 £ E is a reference point. Here and in what follows, the Lebesgue 
measure dx is often omitted. It is convenient for us to define two measures fi 
and v as follows. 

/x(dx) = , , dx, zv(dx) = e '■^Mx. 
a{xj 

The first one has different names: speed, or invariant, or symmetrizable mea- 
sure. The second one is called scale measure. Note that v is infinite iff the 
process is recurrent. By using these measures, the operator L takes a very 
compact form 



d/x du 



i.e., Lf^ae-^if'e^y) (7) 



which goes back to a series of papers by W. Feller, for instance [12J . 

Consider first the special case that M,N < oo. Then the ergodic case 
means that the process has reflection boundaries at —M and N. In analytic 
language, we have Neumann boundaries at —M and N: the eigenfunction g 
of Ai satisfies g'{—M) = g'{N) = 0. Otherwise, in the non-ergodic case, one 
of the boundaries becomes absorbing. In analytic language, we have Dirichlet 
boundary at —M (say): the eigenfunction g of Aq satisfies g{—M) = 0. Let 
us use codes "D" and "N" , respectively, to denote the Dirichlet and Neumann 
boundaries. The corresponding minimal eigenvalues of —L are listed as follows. 

• A^'^: Neumann boundaries at —M and N , 

• A^^: Dirichlet boundaries at —M and A^, 

• A^^: Dirichlet at and Neumann at N , 

• A'^^: Neumann at and Dirichlet at N . 
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We call them the first non-trivial or the principal eigenvalue. In the last 
two cases, setting M = is for convenience in comparison with other results 
to be discussed later but it is not necessary. Certainly, this classification 
is still meaningful if M or is infinite. For instance, in the ergodic case, 
the process will certainly come back from any starting point and so one may 
imagine the boundaries ±00 as reflecting. In other words, the probabilistic 
interpretation remains the same when M, N = 00. However, the analytic 
Neumann condition that lim.x^±oo 9' (x) = for the eigenfunction g of 
may be lost (cf. the first example given in Section [5]). More seriously, the 
spectrum of the operator may be continuous for unbounded intervals. This 
is the reason why we need the L'^-spectral theory. In the Dirichlet case, the 
analytic condition that lima;_!.-|-oo gix) = can be implied by the definition 
given below, once the process goes ±00 exponentially fast. Now, for general 
M,iV ^ 00, let 

D{f)= A^, M,Ar^oo, /Gi2/(-M,iV), 

£/(—M,N) = the set of absolutely continuous functions on (—M,N), 
£/o{-M,N) = {/ G £/{-M,N) : / has a compact support}. 

From now on, the inner product (•, •) and the norm || • || are taken with respect 
to fj, (instead of vr). Then the principal eigenvalues are defined as follows. 



= mf{D{f) 
= inf{D(/) 
A^^ = mi{D{f) 
X™ = inf{D(/) 



/GM)(-M,iV), 11/11 = 1}, (8) 

/ e ^0(0, N), f{N-) = if < 00, 11/11 = 1}, (9) 

/e^(-M,iV), M/) = o, 11/11 = 1}, (10) 

/G^(0,iV), /(0+) = 0, 11/11 = 1}. (11) 



Certainly, the above classification is closely related to the measures /i and u. 
For instance, in the DN- and NN-cases, one requires that /u(0, A^) < 00 and 
^{—M,N) < 00, respectively. Otherwise, one gets a trivial result as can be 
seen from Theorem 12.21 below. 

To state the main result of the paper, we need some assumptions. In 
the NN-case (i.e., the ergodic one), we technically assume that o and b are 
continuous on (— M, N). For A^^ and X^^, we allow the process to be explosive 
since the maximal domain is adopted in definition of A^^ and X^^ . But for 
AND ADD, we are working for the minimal process (using the minimal 

domain) only, assuming that and v are locally finite. 



Theorem 2.2 (Basic estimates [9j) Under the assumptions just mentioned, 
corresponding to each #-case, we have 

{K#y'/A^X* = E^^,^{K*)-\ (12) 
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where 

(kNN)-^ = inf [^,{-M,x)-^ + ^,{y,N)-^]v{x,y)-\ (13) 

x<y '- 

(kDD)"^= inf [K-M,x)"i+z.(y,iV)-i]^(x,y)-\ (14) 

' x<y '- ■' 

sup 1/(0, x) /i(x, A^) (15) 

xe(o, Af) 

= sup /Li(0,x)z^(x,A^). (16) 
xe(o, AT) 

In particular, A* > iff < oo. 

In each case, the principal eigenvalue is controlled from above and below 
by a constant up to a factor 4 which is universal. Among these cases, the 
hardest one is the ergodic case. It may be helpful for the reader to show how 
to write down kJ^^ step by step. 

• We need two parameters, say x and y with x < y. The state space is 
then divided by x and y into three parts: the left-hand part (— M, x), 
the right-hand part {y,N), and the middle one {x,y). 

• Measure the left-hand and the right-hand subintervals by /i and the 
middle one by respectively: 

k = kNN: ^i{-M,x) fi{y,N) u{x,y). 

• Make inverse everywhere: 

: /i(-M,x)-^ Ky,N)-^ i^{x,y)-\ 

• Finally, summing up the first two terms and making infimum with re- 
spect to X < y, we get the answer. 

Every step is quite natural except the second one: why we use /i but not v 
in the first two terms? This is because we are in the ergodic case, is a 
finite measure. If is replaced by v, since z^(— oo,x) and v{y,oo) are infinite 
when M, N = oo, one would get zero for these terms and so the quantity 
is trivial. A sensitive point here is that we use plus, rather than maximum 
in the last step. Otherwise, even though the resulting bounds are equivalent 
to ours but it then would produce a factor 8 rather than 4 as we expected. 
We have thus completed the first, the most important quantity k^'^. To get 
K^^, simply apply the rule: exchanging the codes D and N simultaneously 
in k"^ leads to the exchange of the measures /x and u in the formula. Let 
us now examine (|14p more carefully. When N = oo and iy{y, oo) = oo, the 
second term in the sum of (fT4|l disappeared. In other words, the boundary 
condition D on the right endpoint is replaced by N. Then the variable y is 
free and so can be removed. Therefore we obtain formula (11511. We remark 
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however that the relation between A and k remains the same even if 
v{y^co) < oo. From (fT5|) . using again our rule, we obtain ([T6|). We mention 
that ()16p can be formally obtained from ()13p by removing the second term in 
the sum. Actually, (jl6p is formally a reverse of (jlSp . and so is somehow an 
easy consequence of (fT5]) . 

2.3 Short review on the known results 

It is the position to say a little about the history of the topic. Clearly, we are in 
the typical situation of the Sturm-Liouville eigenvalue problem (1836-1837). 
From which, we learn the general properties of the eigenfunction: the existence 
and uniqueness, the zeros of the eigenfunction, and so on. Except some very 
specific cases, the problem is usually not solvable analytically. This leads to 
the theory of special functions used widely in sciences. The estimation of the 
principal eigenvalues is usually not included in the Sturm-Liouville theory but 
is studied in harmonic analysis (especially for A^^). To see this, rewrite (llip 
as the Poincare inequality 

A°^||/f ^ Z?(/), /(0) = 0. 
More general, we have Hardy's inequality 

pN 

\\f\\U,)^Ap l/Te^' /(0)=0, p>l 

J-M 

where Ap denotes the optimal constant in the inequality. Certainly, A2 = 
^^DN^ i_ This was initialed, for the specific operator L = x^d^/dx^, by G.H. 
Hardy |16j in 1920, motivated from a theorem of Hilbert on double series. To 
which, several famous mathematicians (H. Weyl, F.W. Wiener, I. Schur, et 
al.) were involved. After a half-century, the basic estimates in the DN-case 
were finally obtained by several mathematicians, for instance B. Muckenhoupt 
(1972). The reason should be now clear why (llSp can be so famous in the 
history. The estimate of A^^ was given in Maz'ya (1985). In the DD-case, the 
problem was begun by P. Gurka [H] around 1989 and then improved in the 
book by Opic and Kufner (1990) with a factor ~ 22. In terms of a splitting 
technique, the NN-case can be reduced to the Muckenhoupt 's estimate with 
a factor 8, as shown by Miclo (1999) in the context of birth-death processes. 
A better estimate can be done in terms of variational formulas given in [H 
Theorem 3.3]. It is surprising that in the more complicated DD- and NN-cases, 
by adding one more parameter only, we can still obtain a compact expression 
(I13p and (I14p . Note that these two formulas have the following advantage: the 
left- and the right-hand parts are symmetric; the cases having finite or infinite 
intervals are unified together without using the splitting technique. 
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2.4 Motivation and application 

Here is a quick overview of our motivation and application of the study on this 
topic. Consider the ip^-model on the d-dimensional lattice Z'^. At each site i, 
there is a one-dimensional diffusion with operator Lj = d^/dx^ — u'{xi)d/dxi, 
where Ui{xi) = xf — /3xf having a parameter /? ^ 0. Between the nearest 
neighbors i and j in Z*^, there is an interaction. That is, we have an interaction 
potential H(x) = — ^i^j with parameter J ^ 0. For each finite box A 

(denoted by A d Z'^) and oj G R^'', let H'j( denote the conditional Hamiltonian 
(which acts on those x: Xj. = (jJk for all k ^ A). Then, we have a local operator 

L% = Y,[dii-di{u + H%)di]. 

We proved that the first non-trivial eigenvalue (A, lj) (as well as the log- 
arithmic Sobolev constant o"^(A,a;) which is not touched here) of is ap- 
proximately exp[— /3^/4] — 4dJ uniformly with respect to the boxes A and the 
boundaries lo. The leading rate /3^/4 is exact which is the only one we have 
ever known up to now for a continuous model. 

Theorem 2.3 ([8j) For the ip^-mode\ given above, we have 

inf inf X'^(A,co)^ ini inf cj^(A, w) f«exp r-/3V4-clog /3l - 4dJ, 

where c G [1,2]. See Figure 1. 




= 2dJ 



The figure says that in the gray region, the system has a positive principal 
eigenvalue and so is ergodic; but in the region which is a little away above the 
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curve, the eigenvalue vanishes. The picture exhibits a phase transition. The 
key to prove Theorem 12.21 is a deep understanding about the one-dimensional 
case. Having one-dimensional result at hand, as far as we know, there are at 
least three different ways to go to the higher or even infinite dimensions: the 
conditional technique used in ; the coupling method explained in [H] Chapter 
2]; and some suitable comparison which is often used in studying the stability 
rate of interacting particle systems. This explains our original motivation and 
shows the value of a sharp estimate for the leading eigenvalue in dimension 
one. The application of the present result to this model should be clear now. 



3 Sketch of the proof 

The hardest part of Theorem 12.21 is the assertion for . Here we sketch its 
proof. Meanwhile, the proof for is also sketched. The proof for the first 
assertion consists mainly of three steps by using three methods: the coupling 
method, the dual method, and the capacitary method. 



3.1 Coupling method 

The next result was proved by using the coupling technique. 



Theorem 3.1 (Chen and Wang (1997)) For the operator L on (0, cxd) with 
reflection at 0, we have 



Ai = A^^ ^ sup inf 



b' 



af" + {a' + h)f 
f 



(17) 
(18) 



^ = {/ G <r2(0,oo) : /(O) = 0, /|(o,oo) > 0}. 

Actually, the equality sign holds once the eigenfunction ot Ai belongs to C^. 

We now rewrite the above formula in terms of an operator, Schrodinger 
operator L5. 



Ai = sup inf 



b' 



af" + {a' + h)f 



Lsf 



f 



f 

(x) =: Xs, 



(x) 



= sup inf 

d^ d 
= a(x)^ + (a'(x) + 6(x)) — + 6'(x). 



dx 



(19) 
(20) 
(21) 



The original condition vr(/) = in the definition of A^^ means that / has 
to change its sign. Note that / is regarded as a mimic of the eigenfunction 
g. The difficulty is that we do not know where the zero-point of g is located, 
hi the new formula ()20p . the zero-point of / S ^ is fixed at the boundary 0, 
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the function is positive inside of the interval. This is the advantage of formula 
(|20p . Now, a new problem appears: there is an additional potential term b'{x). 
Since b'{x) can be positive, the operator Lg is Schrodinger but may not be an 
elliptic operator with killing. Up to now, we are still unable to handle with 
general Schrodinger operator (even with killing one), but at the moment, the 
potential term is very specific so it gives a hope to go further. 



3.2 Dual method 

To overcome the difficulty just mentioned, the idea is a use of duality. The 
dual now we adopted is very simple: just an exchange of the two measures /i 
and u. Recall that the original operator is L — by ([7]). Hence the dual 

operator takes the following form 

L* = — — = —— , (22) 
d/x* di^* dz^ d/i' 

d^ d 
L* = a{x)—j + (a'{x)-b{x)) — , xG(0,oo). (23) 
dx^ ' dx 

This dual goes back to Siegmund (1976) and Cox & Rosier (1983) (in which 
the probabilistic meaning of this duality was explained), as an analog of the 
duality for birth-death process (cf. [9] for more details and original references). 
It is now a simple matter to check that the dual operator is a similar transform 
of the Schrodinger one 

L* = e^Lse-^. (24) 



This implies that 



Lsf L*r 



f f* ' 

where /* := e*^/ is one-to-one from ^ into itself. Therefore, we have 
Xs = sup inf ^ {x) = sup inf ^ J (x) = A*°°, 

where the last equality is the so-called Barta's equality. 
we have thus obtained the following identity. 



Proposition 3.2 Ai = As = A*°°. 

Actually, we have a more general conclusion that Ls and L* are isospectral 
from L^(e*^dx) to L^(e~*^d2;). This is because of 

/ e^fLs9 = I e-^ie^f){e^Lse-^)ie^'g) = J e-^fL^g*, 

and Ls and L* have a common core. But L on L'^{fJ-) and its dual L* on 
-L^(e~*^dx) are clearly not isospectral. 
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The rule mentioned in the remark after Theorem 12.21 and used to deduced 
(jl4p from (jl3p . comes from this duahty. Nevertheless, it remains to compute 
for the dual operator. 



3.3 Capacitary method 

To compute A^^, we need a general result which comes from a different di- 
rection to generalize the Hardy-type inequalities. In contract to what we have 
talked so far, this time we extend the inequalities to the higher dimensional sit- 
uation. This leads to a use of the capacity since in the higher dimensions, the 
boundary may be very complicated. After a great effort by many mathemati- 
cians (see for instance Maz'ya 1985; Hasson 1979; Vondracek 1996; Fukushima 
& Uemura 2003; and [7]), we have finally the following result. 



Theorem 3.3 For a regular transient Dirichlet form {D,!^{D)) with locally 
compact state space the optimal constant in the Poincare-type in- 

equality 

\\f\^AMD{f), /G^F(^) 

satisfies ^ ^ 4:B^, where || • ||b is the norm in a normed linear space IB 
and 

Bb = sup Cap{Ky llli^lb- 

compact K 

The space IB can be very general, for instance L^^fj,) {p ^ 1) or the Orlicz 
spaces. In the present context, D{f) = J_j^ f'^e^ , &{D) is the closure of 
^^(-M,iV) with respect to the norm || • = +^(/), and 

Cap(ir) = inf {!)(/) : / G <r^(-M,iV), /|x ^ l}. 

Note that we have the universal factor 4 here and the isoperimetric constant 
5b has a very compact form. We now need to compute the capacity only. The 
problem is that the capacity is usually not computable explicitly. For instance, 
at the moment, I do not know how to compute it for Schrodinger operators 
even for the elliptic operators having killings. Very lucky, we are able to 
compute the capacity for the one-dimensional elliptic operators. The result 
has a simple expression: 

Bm= sup [v{-M,xY^ + v{y,NY^] ^\\l(^^^y)\\n. 

-M<x<y<N 

It looks strange to have double inverse here. So, making inverse in both sides, 
we get 
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Applying this result to B = L^{i^), we obtain the solution to the DD-case: 
ADD ^ ^-1 



-M<x<y<N 



3.4 The final step 



Applying the last result to the dual process and using Proposition 13.21 we have 
not only 

(^*DD)~V4 ^ = Xs = A*D° ^ 

but also 

' x<y '- 

= inf [^(-M,x)-i+/i(y,A^)-^] y{x,y)-^ 

x<y '- -' 

This finishes the proof of the main assertion of Theorem 12. 2[ 
3.5 Summery of the proof 

Here is the summery of our proof. First, by a change of the topology, we reduce 
the study on Emax to A^^. Then, by coupling, we reduce to Xs- Next, by 
duality, we reduce to A*dd. "We use capacitary method to compute A*dd_ 
Finally, we use duality again to come back to X^^. Recall that our original 
purpose is using Ai = X^^ to study the phase transition, a basic topic in the 
study on interacting particle systems (abbrev. IPS). It is very interesting that 
we now have an opposite interaction. We use the main tools (coupling and 
duality) developed in the study on IPS to investigate a very classical problem 
and produce an interesting result. 

4 Improvements 

The basic estimates given in Theorem 12.21 can be further improved. For half- 
line at least, we have actually an approximating procedure for each of the 
principal eigenvalues. Refer to (BJ [9] and references therein. Moreover, one 
may approach the whole line by half- lines. Here we consider an additional 
method but concentrate on A^^ and A^^ only. As will be seen soon, the 
resulting bounds are much more complicated, less simple and less symmetry, 
than those given in Theorem 12. 2i 

Let us begin with a simper but effective result. 



Proposition 4.1 We have 

ADD ^ (-DD)-l ^ (^DD)-l 
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and 



where 

^ x<y \ J 



;NN^-l 



1 



v{-M,x) 

-1 I r Ar\-1 



x\L'{x,y)+ i^{dz) 



M 



1 



K-M,x) 



H 2 ^7V 

X 
2 

+ / i/(dz) 
Jy 



u{y,N) 



n 2n -1 



2n -1 



Note that if v^—M, N) < oo which is not assumed in Proposition 14. 11 then 
the last two terms in {• • • } in the expression of (k^^) can be written as 



ui-M,x)-^ ^x{dz)u{-M,zY+u{y,N)-^ fi{dz)v{z,NY . 



M 



N 



Otherwise, this expression may be meaningless. Similar comment is meaning- 
ful for 

Proof. Fix X < y. Applying A^^ ^ D{f)/ fi[f'^^ to the test function 

Hy,N) 



f{z)={ u{-M,x) 
iy{z,N), 



v{-M, z A x), z ^y 
z^y. 



we obtain A^^ ^ (k^^) ^ . By duality, we obtain the assertion for . Refer 
to the remark after the proof of [H Theorem 8.2] for more details. □ 

To improve the lower estimate Theorem 12.21 we need more work. For a 
given / G '^{-M,N) with f\(-M,N) > 0, define 

e-^(-)dx / ^,z^e, (25) 

-M J X ^ 

h+{z) = h+{z) = Jfi{l^e,.)f)Hz,N))= e-^(^)dx / z>e, (26) 



i.e. (by exchanging the order of the integrals), 

h-iz) = fi{fi^{-M,-Az))=fi(^fi^{-M,-)l^^M,z))+f^{n(z,e))H-M,z), z ^ 9, 
h+{z)=fi{fui-\/z,N)) =/i(/i/(-,Ar)l(,,^)) + /x(/l(e,,))i/(z,iV), z>9, 
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where x Ay = min{x, y}, xV y = max{x, y}, and 9 = 9{f) G (— M, A^) is the 
unique root of the equation: 

h-{e) = h+{9) 

provided < oo. Next, define 11^ {f) = / f ■ 



Theorem 4.2 (Variational formula) Let a and b be continuous and a > 

on {-M,N). 

(1) Assume that u{—M,N) < oo. Using the notation above, we have 

A°°=sup[ inf //-(/)(z)-ij/\[ inf //+(/)(z)-ij, (27) 

where = {/ G '^(-M,iV) : / > on (-M,iV)}. 

(2) Assume that fi{-M,N) < oo. Then ([27]) holds replacing by A^^ 
provided in definition of /i^, p, and are exchanged. 



Proof. By duality, it suffices to prove the first assertion 



a) Without loss of generality, assume that < oo. Otherwise, the 



assertion is trivial. First, we prove Let 

Hz)-- 



h (z), z ^ 9, 
h+{z), z>9, 



Clearly, /i|(_A/,Ar) > and h G 'rf{—M,N) in view of definition of 9. Next, 
note that 



Jx a 

h+'(x) = -e-^(") r — /, h+"{x) = e-^(") 
. f) a 



-C(x) 



a ,}~ a a 



5 rx 



— /-— / 
a la a a 



X <9: 



X > 



Obviously, h'{9 it 0) =0. Since a, b and / are continuous and a > on 
(-M, N), we also have h"{9+0) = h"{9-0) and soke '^'^{-M, N). Therefore, 
by Barta's equality, we have 



A°° = sup inf 

zei-M,N) g 

-Lh 



Lg 



> inf 



ze(-M,Af) h 



(z) 



inf — ; (z' 

ze{-M,e) h~ 



' \ lze(9,N) h+ 
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Now, by ([7]), required assertion follows by a simple computation. 

(b) Next, we show that the equality sign in (|27p holds. The assertion 
becomes trivial if A^^ = 0. Otherwise, the eigenfunction g of A^^ should be 
unimodal (which seems known in the Sturm-Liouville theory and is proved in 
the discrete context [H Proposition 7.14]. Actually, the discrete case is even 
more complex since the eigenfunction can be a simple echelon, not necessarily 
unimodal). By setting f = g and 9 to be the maximum point of g {g'{6) = 0), 
it follows that II^{f)~^ = A^^ and hence the equality sign holds. □ 

We now introduce a typical application of Theorem 14.21 Fix x < y. Define 



/ i^{y,N) 



v{-M, s f\x), s ^y 



s^y 



and set 



inf 

x<y 



sup ii-{r^y){z)]\i\ sup ii+{r^y){: 



By exchanging /u and u, we obtain k^^. Now, by Theorem 14. 2^ we have the 
following result. 



Corollary 4.3 Under assunnptions of Theorem 14.21 we have 



ADD^U^D)-' and ANN^^KNN^-l_ 



We remark that the assumption in part (1) of Theorem 14.21 is necessary 
for DD-case (cf. (|13p ). Recall that (j27p is a complete variational formula 
for the lower estimates of A^^. Starting at /i = / used in Corollary 14. 3^ 
replacing / and h used in Theorem l4.2l bv fn-i and /„, respectively, we obtain 
an approximating procedure from below for A^^. Dually, we can obtain a 
variational formula for the upper estimates of A^^ and an approximating 
procedure from above. Here we omit all of the details. The same remark is 
meaningful for X^'^ , which is especially interesting since here we do not use 
the property that fi{f) = for the test function /. The new difficulty of ([27|) 
is that 9{f) may not be computable analytically. This costs a question to 
prove that kP^ ^ AkP^ which should be true in view of our knowledge on the 
half-line, and is illustrated by examples in the next section. It is noticeable 
that the method works for the whole line and the use of 0(f) is essentially 
different from what used in the splitting technique. Finally, we mention that 
the method used here is meaningful for birth-death processes, refer to [H 
Lemma 7.12]. 
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For convenience in practice, we express used in Corollary 14.31 more 
explicitly. Let i^-{s) = v{—M,s) and v+{s) = u{s,N) for simplicity. Then 



f{s) = rHs) 



and 



h (^) = /^(/i^-l(-M,^)j + fi{fl(z,e)) 



S ^ X 

X ^ s 



z ^ 6 
z > 



(28) 



(29) 
(30) 



We now consider the typical case that Q G \x^y\- Then, 
Hence the equation h~{6) = h^{6) becomes 



i==^(i^+^l(j;,Ar)) +//(^z^+]l(e,y)), Oe [x,y]. 



vMy) 



(31) 



Furthermore, by some computations, we obtain the ratio / f^'^ as follows. 
We have for z: z ^ x ^ 9 ^ y that 



+ V-{z)v^{x) ^J,{x,6), 

and for z: z y 9 that 



(32) 



+ Vi/+(z)zv+(y)/i((9,y). 



(33) 



Note that by (j25|) and (|26|) . /i is increasing on [x, 0] and is decreasing on 
Since /^'^ is a constant on it follows that 



h-{z) h'(9) , /i+(z) /i+(e) 

max — -— = — -— and max — 



ze[x,e] f^'y{z) f^'y{x) 
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By assumption, h {9) = h^{9). Hence 



max ir{f^y)(z) = max II+{f^y)(2 



1 



Thus, for computing k^^, by (f32]) ~ (pl]l . we arrive at 

sup ii~iry)iz)]\j\ sup ii^{ry){z) 



(34) 



z£{-M, 



sup 

ze{-M, x) 



ze{e,N) 



^ ^^{ly^J'^t{-M,x)) + f^[^-^{x, 



ze{y,N) 



(35) 



Finally, let (x* ,y* ,9*) solve equation (j3ip and two more equations modified 
from ignoring its left-hand side and replacing the last two "V" with "=" . 
Then we have 



DD 



1 



^/uIix*) 



^1 ) + f^['^-'^{x*,e*) 



(36) 



5 Examples 



This section illustrates the application of the basic estimates given in Theorem 
and the improvements given in Proposition 14.11 and Corollary 14.31 



Example 5.1 (OU-processes) The state space is M and the operator is 

2 \ dx'^ ax 

This is a typical example of the use of special functions. It has discrete eigen- 
values Xn = n with eigenfunctions (Hermite polynomials) 

5„(x) = (-l)"e-— (e-), n^O. 
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Then, we have (k^^) = Aq = 0, A^^ = Ai = 1 with eigenfunction g{x) = x. 
To compute , noting that the operator, the eigenfunction are ah symmetric 
with respect to and so does k^'^, one can spht the whole line into two parts 
(—DO, 0) and (0, oo) with common Dirichlet boundary at 0. This simplifies the 
computation and leads to (k'^'"^) ^ = (k^'^) ^ « 2.1. Note that g'{x) = 1 but 
limi^Koo {e.^g'){x) = 0. 

For the half-space (0,c«), as we have just mentioned, A^^ = A^^ = 1 
with g{x) = X, (k^^) ^ = (k^^) ^ ~ 2.1. For , the symmetry in the 

whole line is lost. We have A^^ = 2 with g{x) = -1 + 2x'^, {k^^)~^ ^ 4.367 
which is achieved at {x,y) ^ (0.316,1.185). Note that liuix^oog'ix) = oo but 
lima^^oo {e'^g') (x) = 0. 

To study , recall that we can reduce the NN-case to the DD-one 
by an exchange of /i and i/. By Proposition 3.1, we have (rNN)-^ ~ 2.6. 

By Corollary [41 and ([MD, we obtain (k'^^)"^ 1.83 with {x*,y*,9*) ^ 
(0.6405,0.938,0.721194). For the last conclusion, we use a direct search star- 
ting from {x,y) ~ (0.316,1.185) which leads to in the last paragraph. 
The ratio becomes 2.6/1.83 ~ 1.42 < 4. We mention that similar estimates 
can also be obtained by using a different approximating procedure in parallel 
with [21 Theorem 6.3]. Refer to [SJ Footnotes 12 and 14]. 

The following examples are often illustrated in the textbooks on ordinary 
differential equations, see for instance Hartman (1982), §11.1. 



Example 5.2 The equation 

u" + (T^n = (o- / 0) 

has the general solution 

u = ci cos{ax) + C2 sin((Tx). 

From this, it should be clear that for the operator L = /dx"^ with finite state 
space {a, (3), we have 



13 — a) '' \ f3 — a J ^ 

A°-=(;,7^V, .(x)=sm^^(^-") 



2{p-a)J ' ' V2(/3-a) 

A"°=(:,7^V, .(.)=cos^^(^-«) 



2{P-a)J ' ' V2(/3-a) 



18 



Mu-Fa Chen 



The corresponding estimates are as follows. 

2 



/3 — a 



Note that by symmetry, the DD- and NN-cases can be split at 6 = (a+/3)/2 
into the DN- and ND-cases. One can then approach and by using 
the known approximating method for A^'^ and A^^ (cf. [SJ Theorem 1.2]). 
However, as an illustration of Theorem l4.2l and Corollary 14. 31 we now compute 
and kDD. 

Consider first the simpler interval {a, 13) = (0, 1). Since ^ = v = dx, by 
symmetry, one may choose y = 1 — x. Then x < 1/2 and 



(kDD)-'= inf - 
xe(o,i/2) X 



l-2x + x-^ z'^dz + x'^ (l-zfdz 

Jo Jl-x 



-1 



6 



inf , , „ 

xe(o,i/2) 3x - (2x)^ 

32 



(with X = 3/8). 
o 

To compute kP^ , set again y = 1 — x with x G (0, 1/2). Then, the test 
function /^'^ becomes 



V-s A X 



S ^ 1 — X 

s G (1 - X, 1). 



By symmetry again, we have 9 = 1/2. Fix x G (0, 1/2). For convenience, we 
express as (/i, /2): /i(s) = for s G [0, x] and /2(s) = \/x for s G [x, 1/2]. 
Then by (p9]) with z^-(s) = s, we have /i~ = (^h^ (z) , h2 (z)) : 



z e [0, x] 



/ fi{s)sds+ [ f2{s)sds] + Z f ' f2, [x,l/2]. 



' PX 


/■1/2 -1 


. J 2 


/ ^2 



/l2 (z) 

Hence by we have 
h-{z) 



ii~{n{z) 



i 7.3/2 I 1^1/2 



2 15^^' 



/^(^) U(5.(l-.)-x2), 



z G [0,x], 
z G [x, 1/2]. 



Define 



Then 



H{x) 



-ix3/2 + LV2 and 7(z) = F(x) - ^z^. 



H{x) 8 F(x) 8 



< 0. 
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15 ^ 



Hence 7 achieves its maximum at 



Furthermore, 

Note that z*{x) ^ x iff x ^ 5/14. Besides, on the subinterval [x, 1/2], 

h~ (z) / (z) has maximum 1/8 — x^/ 10 by (j34|) . Solving the equation 

f)%x)^/^ = ^-3^-^ (5/14, 1/2), 

we obtain x* ~ 0.436273 and then 

inf sup ^^44 = liz*(x*)) w 0.105967. 

.6(5/14,1/2) ,^1^2 Piz) 



From these facts and (|36p . we conclude that 

~ 1/0.105967 9.43693. 
By the way, we mention that a similar but simpler study shows that 

. r h-{z) 1 

mt sup = -. 

a^e(o,5/i4) ^^1/2 r{z) 8 

This shows that to get a less sharp lower bound 1/8, the computation becomes 
much simpler. It needs to study the extremal case that x = only; the 
corresponding test function becomes = 1. Return to the original interval 
(a,/3), by Proposition 14.11 and Corollary 14.31 we obtain 

8 ^ 9.4369 ^ ^ / y ^ 32 ^2/4 



(/3-a)2 (,S-a)2 \fi-aj 3(/3 - a)2 3V/5-a 

The ratio becomes ^/9.4369 1.13. The same assertion holds if A^^ is 
replaced by A'^^ because of the symmetry. 

It is a good chance to discuss the approximating procedure remarked after 
Corollary 14.31 Here we consider the lower estimate only. Replacing = 
(/i)/2) by one produces a new and then a new 11^ (f) 

which provides a new lower bound. By using this procedure twice with fixed 
6 = 1/2 and x = x* ~ 0.436273, we obtain successively the following lower 
bounds: 

9.80392 9.86193 



20 



Mu-Fa Chen 



Clearly, they are quite close to the exact value of A^'-' and X^^: 

TT^ _ 9.8696 
{P - a)2 ~ (^-a)2- 

Example 5.3 By a substitute u = ze"^^/^, the equation 

u" + bu' + 'yu = (6, 7 are real constants) 

is reduced to 

z" + a^z = (^2 = 7- 62/4). 
From the last example, it follows that the equation has general solutions 

'e-^^/2(ci +C2a;) if 7 = ^2/4 

u = J cie^i^ + C2e^2x if ^ < j,2/4 

e-''^/2 (|ci cos {x^yJ-b'^/A) + C2 sin (x^T-^V^)) if 7 > &V4, 
where ^1 , ^2 are solution to the equation 

f + b^ + j = 0. 

Thus, for the operator L = d2/da:2 + 6d/dx (6 is a constant) with state space 
(0,00), wc have the following principal eigenfunctions 

• g{x) = (2/6 + x)e~^^/2 ^nd g{x) = xe~^^/'^ in ND- and DD-cases, respec- 
tively, when 6 > 0; 

• g{x) = a;e-^^/2 and g{x) = (1 + 6a;/2)e-^^/2 in DN- and NN-cases, 
respectively, when 6 < 0. 

In each of these cases, we have the principal eigenvalue A* = ^2/4 and 
{K*y^ = 52. Moreover, {kP^)'\ {r^^)'^ = b^ /2. Clearly, the lower es- 
timate (k*) ^/4 is sharp in all cases. 



Example 5.4 (Cauchy— Euler equation) Consider the operator 

L = X -p-^ -I- ox—, 
ax^ ax 

where 6 is a constant. By a change of variable x = e^, the equation 

x^u" + bxu' -|- 7U = (^,7 are constants) 
is reduced to the last example: 

d'u ^.du 
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u= < 



Hence the original equation has general solutions 

'x(i-'')/2 (ci + C2logx) if 7 = (1 - 6)2/4 

cia;«i + C2X«2 if 7 < (1 - 6)V4 

2:(i-6)/2 cos (a/7 - (1 - 6)2/4 log x) +C2 sin (^7 - (1 - 6)^/4 log x)^ 

if7>(l-5)V4, 

where ^i, ^2 are solution to the equation + {b — 1)^ + 7 = 0: 

6,6 = (1 - b)/2 ± 7(1-6)2/4-7. 

Here we have used Euler's formula: 

x'^ = e'^ ^''^'^ = cos {^/i\ogx)+i sin ( log x) . 

In particular, 

(1) when 6 = 0, we have solutions 

{^/x{cl + C2logx) if 7 = 1/4 

cix^^ + C2X^'^ if 7 < 1/4 

^/x(^Cl COS — 1/4 loga;)+C2sin (-^7^^174 logx)^ if 7 > 1/4. 

Now, corresponding to 7 = 1/4, we have 

, pvivT 1 , , I \/x if the state space is (0, 00) 

A = -, g[x) = < 

4 log ^/x if the state space is (1, 00). 

The first case is the original Hardy's inequality. Corresponding to 7 = 1/4 
again but for state space (l,oo), we have 

A^^ = ^, g{x) = V^{logV^-l). 

Here lim.^oo {e^g'){x) = lim^^oo ^'(x) = 0. We have (k^^)"' = 1, 

^, (r^^) ^ = 1/2, respectively. The lower estimate ^/4 is sharp 
in each case. The DN-case is actually a special one of the last example. 

(2) When 6 = 1, for finite state space (1, N) with Dirichlet boundaries, we 
have 

In particular, 

,DD TT \- _ . / TT 
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Next, for Neumann boundaries, we have 



2 



In both cases, we have , (k^N) = (4/logA^)^ Besides, we have 



,DN ' ^ ^^ ' ^ . / vr 



In these cases, we have (k^^) ^, (k^^) ^ = (2/logA^)^. Note that the present 
case can be reduced to Example 15.21 under the change of variable x = e^, the 
results here can be obtained from Example 15.21 replacing (a — by log^ A^. 
In view of this, we also have 



^ ^ ^ ^ Slog^iV ^ ^ log^iV 



6 Appendix 

The next result is a generalization of [9l Proposition 1.2]. 



Proposition 6.1 Let Pt{x,-) be symmetric and have density pt{x,y) with re- 

1/2 

spect to fj,. Suppose that the diagonal elements ps{-, •) G L-^^^ {jj.) for some s > 
and a set of bounded functions with compact support is dense in Then 

■^0 — £^max- 



Proof. The proof is similar to the ergodic case (cf. [U Section 8.3] and [H 
proof of Theorem 7.4]), and is included here for completeness. 

(a) Certainly, the inner product and norm here are taken with respect to 
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//. First, we have 

Ptix,K) = PsPt-slKix) 



dP (x ■) 

fi{dy) 7 ' ^ {y)Pt-stK{y) (since < fi) 



dPsjx,-) 
d/i 



Pt-slK 



dPs(x, 

fi{ IxPt^s — 1 (by symmetry of Pt) 



dfi 



P- 



dPs{x, 



t-s- 



dfj, 



dPs{x,- 



d/i 



(by Cauchy-Schwarz inequality) 
-^o(t-s) |-|-jy ^2_g^pQj^gj^|.-j^j convergence) 



lx{K)p2s{x,x) 



(by [S (8.3)]). 



By assumption, the coefficient on the right-hand side is locally /x-integrable. 
This proves that Smax ^ Aq- 

(b) Next, for each / E with ||/|| = 1, we have 



\Ptff = if, P2tf) (by symmetry of Pt] 



supp (/) 



f,idx)P2t\f\ix) 



supp (/) 



/x(dx)P2t(a;,supp(/)) 



supp (/) 



/i(dx)c(j:,supp(/))e^ 



ax*- 



=: Cfe 

The technique used here goes back to Hwang et al. (2005). 

(c) The constant Cj in the last line can be removed. Following Lemma 
2.2 in Wang (2002), by the spectral representation theorem and the fact that 
= 1, we have 



;>oo 

\Ptff= / e-'^'d{ExfJ) 
Jo 

/ e-^^'d{Exf,f) 
Jo 



t/s 



(by Jensen's inequality) 



t > s. 



Note that here the semigroup is allowed to be sub-Mar kovian. Combining this 

1 2 ^ -2em, 



with (b), we have ||P./f ^ Cfe 



Letting f — t- oo, we obtain 
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first for all / G .J€ and then for all / € L'^ifJ-) with ||/|| = 1, because of the 
denseness of in L'^{fi). Therefore, Aq ^ Emax- Combining this with (a), we 
complete the proof. □ 

The main result (Theorem [22]) of this paper is presented in the last section 
(section 10) of the paper [9], as an analog of birth-death processes. Paper [9], 
as well as |8] for (/9^-model, is available on arXiv.org. 
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